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Potentiality of the presence of traversable wormholes in the outer/inner regions of the halos of 
galaxies, situated on the Navarro-Frenk-White (NFW) density profile and Universal Rotation Curve 
(URC) dark matter models have been investigated recently [27H^ [49H5^ . Since this covers our 
own galaxy also as a possible home for traversable wormholes it prompts us to further the subject 
by considering alternative density distributions. From this token herein we make use of the Einasto 
model [TJl [13| [35] to describe the density profiles for the same purpose. Our choice for the latter 
is based on the fact that theoretical dark matter halos produced in computer simulations are best 
described by such a profile. For technical reasons we trim the number of parameters in the Einasto 
profile to a possible minimum. Based on such a model it is shown that traversable wormholes in the 
outer regions of spiral galaxies are possible while the inner part regions prohibit such formations. 
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I. INTRODUCTION 

Although Flamm’s work on the wormhole physics dates back to 1916, in connection with the newly found 
Schwarzschild solution [18], wormhole solutions were firstly considered from physics standpoint by Einstein and 
Rosen (ER) in 1935, which is known today as Einstein-Rosen bridge connecting two identical space-times [14]. Then 
in 1955 Wheeler used ”geons” (self-gravitating bundles of electromagnetic fields) by giving the first diagram of a 
doubly-connected-space [58]. Wheeler added the term ’’wormhole” to the physics literature, however he defined it at 
the quantum scale. After that, first traversable wormhole was proposed by Morris-Thorne in 1988 [36]. Then Morris, 
Thorne, and Yurtsever investigated the requirements of the energy condition for wormholes m- After while, Visser 
constructed a technical way to make thin-shell wormholes which thoroughly surveyed the research landscape as of 
1988 [471I6TH63] . After this, there are many papers written to support this idea [22l(331 [341(43] . 

In 1933 E. Zwicky, was the first astronomer to propose the existence of the dark matter (DM) after that he observed 
the motions of the galaxies in the Coma Cluster, a galaxy cluster roughly 320 million light-years away and nearly 2 
light-years across, and found that it moved too rapidly [65]. He noted that its speed of revolution, which depends 
on the weight and position of the objects inside, implied the cluster had much more mass than he could see. He 
decided there must be a hidden ingredient, known as dark matter, that caused the motions of these galaxies to be 
so large. In 1978, American astronomer Vera Rubin looked at individual galaxies and she realized thet the outer 
stars were spinning as quickly as the interior stars, and sometimes faster [59]. It is concluded that the galaxies 
must be surrounded by a halo of matter which is not seen. Theorists have postulated many candidates for what 
DM might be. But the most widely accepted hypothesis is the weakly interacting massive particles (WIMPs), a 
particle that interacts with gravity but not light, hence its invisibility to us [7]. Neither WIMPs nor any other particle 
that could successfully explain DM exist in the standard model of particle physics, the theoretical framework that 
has successfully predicted nearly all the phenomena in the universe. Thus, the identity of the DM particle remains 
one of the outstanding mysteries in particle physics and cosmology. According to the latest results from the Planck 
satellite [46], a mere 26.8% DM percent of the universe which plays a central role in modeling of cosmic structure, 
galaxy formation and evolution and on explanations of the anisotropies observed in the cosmic microwave background 
(CMB), 4.9% is made of ordinary matter which is matter composed of atoms or their constituents and 68.3% is dark 
energy (DE) that tends to accelerate the expansion of the universe [8] (20] [21] . Certain types of DM can actually help 
the formation of wormholes. Eor example, dissipative DM in which DM particles forming galactic halos are presumed 
to have dissipative self-interactions can aid the formation of black holes by dissipative collapse, and subsequent to 
wormholes. Now, in this types of models one can consider either DM collapsed to a disk, or a quasi-isothermal profile 
[I6l[l7l[l9l[26l[38]. 
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Recently, Rahaman et al. have used both the NFW mi EH EH 


p{r) = 



and URC DM density profiles [4] 


p{r) 


Porp 

(r + ro)(r2 + rg) 


( 1 ) 
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where pQ, vq and are all constants, to show the potentiality of the presence of traversable wormholes in the 
outer and the inner part, respectively, of the halo. 

In this paper, following similar line of thought, we study the Einasto model nainiES] which is the special function 
that arranges the finest total fit to the halo density profiles. It tells how the density p of a spherical stellar system 
changes with distance r from its center |23]: This is expressed by the density profile 


p(r) oc 




( 3 ) 


with constants A and a. 

The bounded parameter a which is called Einasto index, adjusts the degree of curvature of the profile [9[ [53] 



FIG. 1: Einasto profile in log-log plot (the parameter a 


eontrols the degree of eurvature of the profile) from Eqn. . 


p(r) = po exp 


a ro 
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simply named as the Einasto model in the present paper which is plotted in Eig. Eor the Milky Way galaxy, we 
have ro = 9.11 kpc and po = 5 x 10“^^(ro/8.6 kpc)“^gcm“^ j4l[32]. 

The slopes of the NEW profile has the range from the inner/outer ( r ^ 0 / r ^ oc ) are depicted in Eigs. 
and below. On the other hand, the values are different for the Einasto profile which are included as limiting cases 
of isothermal and Gaussian, (a = 0, a = 2, respectively). The difference between the Einasto and NEW profiles 
appears at small and large radial distances. After all, on scales of concern for gas or stellar dynamics and strong 
gravitational lensing the Einasto profile with a = 0.25 is very akin to the NEW profile, and on scales of concern for 
weak gravitational lensing the Einasto profile with a = 0.2 is quite similar to the NEW profile. Galaxy mass halos 
typically have a = 0.16 as referred in Eigs. [^and §m- 
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Einasto’s model has been used in the description of many systems, such as galaxies, and halos. It is used to show 
the DM profiles of simulated halos as well as the NEW profile [UliES]. The Einasto profile has a finite (zero) central 
slope, in contrast to the NEW profile which has a divergent (infinite) central density. Due to the restricted resolution 
of N-body simulations, it is unclear yet which model description fit well to the central densities of simulated DM 
halo. The observations of both the Milky Way and some galaxies may be compatible with the NEW profile for the 
DM halo [24l [3TJ HOI |57l ED] . There are also self-interacting dark matter (SIDM) simulations and baryon feedback 
simulations [illlllSl]. The SIDM is attractive because it offers a means to lower the central densities of galaxies 
without destroying the successes of CDM on large scales. Cosmological simulations that contain only CDM indicate 
that DM haloes should be cuspy and with (high) concentrations that correlate with the collapse time of the halo. One 
possible answer is feedback in which the expulsion of gas from galaxies can result in lower DM densities compared 
to dissipationless simulations, and thus bring CDM models in line with observations [54lEl]- Eurthermore baryon 
feedback has most significant effects on the inner density profiles which are typically shallower than isothermal and 
the halo concentrations tend to be lower than in the absence of baryons. Eor instance the strong feedback models 
reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback [6]. 

The DM profile of smaller galaxies tend to have flatter distributions of DM in the central region, known as the 
cuspy halo problem mills]. Our purpose in this paper is to search for possible traversable wormholes in a spiral 
galaxy specified by the Einasto density profile. Wormholes are solutions to Einstein’s equations that connect different 
universes or tunnels through different parts of the same universe. Our primary aim is to consider the Milky Way 
galaxy as a test bed for this venture. To this end we choose the parameter vq of Milky Way to coincide with the throat 
of the traversable wormhole whereas the second parameter a is in a restricted range. Erom the Einstein equations we 
derive the equations for the density and pressures in terms of the metric functions and their derivatives. Although 
our analytical functions are expressed in terms of the Whittaker functions [T] the detailed analysis are not imperative. 
It suffices for us to check the null-energy condition and satisfaction of the flare-out conditions [36j [37] . The same 
technical handicap prevents us to investigate the stability of the resulting traversable wormhole that may exist in 
the Milky Way galaxy. It is observed that at the central region Einasto profile does not give a traversable wormhole 
solution because violation of null energy condition is not satisfied there. But in the outward region we obtain possible 
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FIG. 3: Comparison of NFW (solid black lines) and Einasto (colored lines) profiles. In this example the halo eoneentration C 200 = 7.1 

m 



FIG. 4-- Rotational veloeity curve for Milkway (the eontribution to the rotation curve as predicted from different models for the stellar 
bulge (blue), stellar disk (green) and gas (blaek). We assume a distanee to the galaetie eenter ro = 8 kpe \24V 


traversable wormholes. 

The paper is organized as follows. In Sec. II we present the Einstein equations for wormhole space-time and we 
find their solutions under the Einasto DM profile. In Section III we discuss our results and conclude the paper. 


II. TRAVERSABLE WORMHOLES UNDER THE EINASTO DM PROFILE 

The spherically symmetric and static wormhole geometry which was proposed by Morris-Thorne is given by the 
line element [36] 


1 - 


b{r) 


ds^ = -e^rOdP + 


dr^ _l_ r‘^(d0‘^ + sin^ 0 d(jp‘) 


( 5 ) 
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FIG. 5: Rotational curves v(r) for some galaxies from THINGS survey together with fits. Blue thiek curves show best fits by DM 
distributions of semi-degenerate configurations, while magenta thin curves show Einasto profile best fits m 


where the function 6(r) defines the spatial shape function, and /(r) stands for the redshift function. The radial 
coordinate r ranges from +oc to 6(ro) = ro(minimum radius condition that any wormhole must satisfy). The minimum 
value of the r is ro, called the throat of the wormhole. The proper radial distance is given by 



Properties of such traversable wormholes are listed below [36l |37] ; 

• Spherically symmetric, static metric, 

• Satisfy the Einstein field equations, 

• No event horizon, i. e. ^ 0, 

• ’’Throat” r = vq connects two asymptotically flat space-time regions, (for Milky Way galaxy ro = 9.11 kpc) 

• The stress-energy tensor (violating the null-energy condition (NEC) with p Pr < 0, where p is the energy 
density and Pr the radial pressure), 

• Physically reasonable construction materials, (for Milky Way galaxy po = 5 x 10“^^(ro/8.6 kpc)“^ g cm 

• Bearable tidal gravitational forces, 

• Satisfy the flare-out condition at the throat {b'{rth) < 1): while b{r) < r near the throat, 

• Einite and reasonable crossing time, 

• Stable against perturbations. 

The Einstein field equation are given by 

= 87rT/f 


( 7 ) 
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where Tl^ is the stress-energy tensor, is the Einstein tensor with the units in which c = G = 1. Field equations 
relate space-time curvature to matter and energy distribution. The non-zero components are 


G 



( 8 ) 



8(1-‘)£ 

r r 


( 9 ) 




f" + +f 




1\ f b'r-b y 
r / \ 2r(r — b) j ’ 


( 10 ) 


in which a prime is 

It is assumed that DM is expressed in the form of general anisotropic energy-momentum tensor 13110], as follows 

= (p-hPr)u^U^ +Pr9y + {Pt -Pr)V^Vu^ (H) 

in which = —\r]^r]^ = —1, and where pt is the transverse pressure, Pr is the radial pressure and p is the energy 
density. A possible set of and are given by 0, 0, 0) and = (0, 0, ^ e )• once 

we restrict ourselves entirely to a collisionless cold DM the pressure terms should not be taken into account. Our 
choice implies accordingly that we aim a more general anisotropic fluid, at least for our future interest. The only 
non-zero components of are 


p = -p, 

(12) 

P = Pr, 

(13) 

II 

II 

(14) 


Furthermore we need to impose some constraints in order to solve Einstein field equations. One of them is the 
tangential velocity, (for a justification of this condition we refer to the book by Chandrasekhar [5l [3Q]i 


= 


(15) 


which is responsible to fit the flat rotational curve for the DM. Other one is proposed by Rahaman et. al m that 
the observed rotational curve profile in the DM region is given by 


= ar exp(—/cir) -h /3[1 — exp(—/c 2 r)] 


(16) 


in which a, /3, /ci, and k 2 are constant positive parameters. It is illustrated in Figj^ With this choice it is guaranteed 
that for r —oo we obtain = [3 =constant which is required for flat rotation curves. 

From the Eqns.(15) and (16), one obtains the redshift function as follows 


f{r) = - 


O 

a r 


a 


2aP 


2fcie(2feir) 4^2g(2feir) ki + k2 

+/32 ln(r) + 2(i‘^Ei{l, fear) - 2k2r) + D. (17) 

with the exponential integral Ej [25| and integration constant D. For large r, it is required also that = 

By using the Eqns. (15) and (16) and the Einasto density profile (Eqn. 0), the Einstein field equations give us 
the tedious form of the shape function 


b{r) = ?T_d5^exp(-) 
a a 


(r/ro)^r“ 


-a \ -3/a 


3(cr “h 3)(2cr T 3) 


Q,3g32(3/a-l-(a+3)/2a) f (f/roYr 


-a\ 3/a 


(18) 
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FIG. 6: Proposed rotational velocity with values of the parameters as ki = k 2 = !;« = 0.0006 and /3 = 0.00082 m 



FIG. 7; f(r) versus r (ki = k 2 = ^,ol = 0.0006 and (5 = 0.00082 ]50^ ) 
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a J a 

where the integration constant is chosen as (in order that 6(ro) = tq is satisfied) 

-ro 


a 


Ci = 


6(^2 + 21a + 27 


(yl7rr^exp(l/a)2(3«-3)/2«)Q,i+(«+3)/2a^ 


( 19 ) 


[(a(5 + a)WhittakerM{—{a + 3)/2(a, (2(a + 3)/2(a, 2/(a) 


+ (3 + a)WhittakerM{{a + 3)/2a, {2a + 3)/2a, 2/a) — 6a — 9]. 

Note that in this expression WhittakerM stands for the Whittaker function [1]. The second order differential 
equation satisfied by b{r) in our present analysis is given by 


6" = ?^ 
r 


1 - ( - 

^0 


( 20 ) 


which follows from Eq.(12). The redshift function /(r) versus radial coordinate r is plotted in Figj^ 



FIG. 8: The shape function b versus radial coordinate r with the parameters ro = 1 and po = 0.0001 


It can be checked that the wormhole throat condition of b{ro) = ro is satisfied. It is shown in Figs. 10 and 11 For 
the Milky Way galaxy, we have ro = 9.11 kpc and po = b x 10“^^(ro/8.6 kpc)“^ g cm 


-3 


in order that ro 

is greater than the horizon of the supermassive black hole. This is necessary for the wormhole to be traversable. The 
choice of ro as throat radius can be justified from the fact that ro plays the role of the core radius for the galaxy. For 
this reason we rely on Fq.Q since for r > tq there is a turning point for the density distribution. We admit, however, 
that the exact choice of ro as the throat radius, unless a better choice is made, is open to criticism. Furthermore, the 
flare-out condition {b' < 1) 


b' = Snpr^ = Sirr^po exp 


a ro 


( 21 ) 
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FIG. 9: The radius of the throat is obtained where b(r)-r euts r axis or the values ro = 1 and po = 0.0001 



FIG. 10: b' versus r for Milky Way (ro = 9.11 kpe and po = 5 x 10 ‘^^{ro/S.6 kpc) ^ g cm ^ ^ \32\ \4^ ) 


is also satisfied for the Milky Way galaxy 6'(ro) = 0.98 x 10“^^ < 1 (Figp^ and also all spherical stellar systems 
as shown in Fig. In addition, both the star-count and kinematic data of the Milky Way stellar halo are well- 
represented by an Einasto profile with index a ~ 0.5 and effective radius ~ 20 kpc, if the dark halo has a flat rotation 

curve (Figs, [^and Em ■ _ 

Finally, for the sake of completeness the radial and lateral pressures are calculated easily from Eqs.(8p4) as 
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FIG. 11: h' versus r (ro = 1 and po = O.OOOll 
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(22) 

r r J 

W / h'r -h 


(23) 


in which /(r) and b{r) are to be substituted from the solution (Eqns. 17- p!9). 

Note that the null energy condition (p+Pr < 0) is violated so that the flare-out condition for a wormhole is satisfied 
according to Fig. Let us note also that those pressures have not been used in the present study but it can be 
anticipated that for a future stability analysis they will be needed. 


III. CONCLUSION 

Existence of traversable wormholes is an important problem in physics both at micro and macro scales. There is no 
doubt that together with black holes wormholes constitute the most interesting objects in our entire universe. On the 
one hand collision of ultra-energetic particles may create these objects while at the other extreme they form naturally 
in galactic systems of cosmic dynamics. Once formed, in order to survive traversable wormholes must satisfy certain 
criteria of stability, otherwise they will decay instantly into more stable objects. Besides stability energy conditions 
are also of prime importance concerning traversable wormholes. This demands negative energy density which is not 
available in classical physics although there are rooms for it in the quantum domain. Let us add also that a recent trend 
is to modify the geometry of throats and obtain thin-shell wormholes [22j [43] supported by positive total energy, for 
this see [33l [34] and references cited therein. However, the minimum requirement of violationg null-energy condition 
(NEC) provides that p + < 0, without specifying the (-) condition for the energy density. As a matter of fact NEC 

provides us a jumping board into the realm of traversable wormholes. Stability, on the other hand, admittedly and for 
completely technical reasons remains open in our present study. This item may be considered as a separate subject 
in our future analysis. 

This work is motivated mainly by Ref. [50|, which discusses the possible formation/existence of traversable worm- 
holes in our own galaxy. Milky Way. Their conclusion that traversable wormholes may exist both at the inner and outer 
parts of our spiral galaxy has been revised by using a different density distribution. Namely, we adapt the Einasto 
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FIG. 12: Null energy eondition (p-\-Pr < 0) for the values of parameters ro = 1 and po = 0.0001 


density profile [T2j [13] which differs much from the Navarro- Frenk- White (NFW) profile employed in [4Tl|50]. Our 
analysis suggests accordingly that at the central region since the NEC is not violated, we do not get any traversable 
wormholes. This explains also the choice of the throat outside of the central region. Yet at the outer regions, the NEC 
is violated so that the fiare-out condition for traversable wormhole is satisfied and we expect formation of traversable 
wormholes. This result is not specific to Milky Way alone but is valid for any spiral galaxy once it obeys a mass 
distribution given by Einosto profile. It is worthwhile therefore to check other mass distributions and see the resulting 
consequences. In conclusion we wish to comment that the DM lurking in the innermost part of the Milky Way [24|, 
may support the existence of traversable wormholes that contains our solar system as well. 
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